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ABSTRACT 

In  this  paper,  we  investigate  the  dynamics  of  a  nonlinear  delay  differential  equation  model  for  passive  mode-locking  in 
semiconductor  lasers,  when  the  delay  model  is  seeded  with  parameters  extracted  from  the  gain  and  loss  spectra  of  a 
quantum  dot  laser.  The  approach  used  relies  on  narrowing  the  parameter  space  of  the  model  by  constraining  the  values  of 
most  of  the  model  parameters  to  values  extracted  from  gain  and  loss  measurements  at  threshold.  The  impact  of  the  free 
parameters,  namely,  the  linewidth  enhancement  factors  that  are  not  available  from  the  gain  and  loss  measurements,  on 
the  device  output  is  then  analyzed  using  the  results  of  direct  integration  of  the  delay  model.  In  addition  to  predicting 
experimentally  observed  trends  such  as  pulse  trimming  with  applied  absorber  bias,  the  simulation  results  offer  insight 
into  the  range  of  values  of  the  linewidth  enhancement  factors  in  the  gain  and  absorber  sections  permissible  for  stable 
mode-locking  near  threshold.  Further,  the  simulations  show  that  this  range  of  permissible  values  progressively  decreases 
with  increasing  bias  voltage  on  the  absorber  section.  This  is  important  for  telecomm  and  datacom  applications  where 
such  devices  are  sought  as  pulsed  sources,  as  well  as  in  military  RF  photonic  applications,  where  mode-locked  diode 
lasers  are  used  as  low  noise  clocks  for  sampling. 

Keywords:  Semiconductor  lasers,  quantum  dots,  delay  differential  equation,  passive  mode-locking,  linewidth 
enhancement  factor. 


1.  INTRODUCTION 

Quantum  dot  mode-locked  lasers  (QDMLLs)  are  attractive  optical  sources  for  applications  requiring  pulse  trains  with 
low  timing  jitter  and  low  energy  per  bit,  such  as  optical  interconnects  in  multi-processor  architectures1’ 2.  In  order  to 
optimize  QDMLLs  for  these  precision  applications,  it  is  necessary  to  accurately  predict  device  operation  using  realistic 
parameters  characteristic  of  the  underlying  semiconductor  gain  material,  under  realistic  operating  conditions.  To  this 
end,  certain  properties  distinctive  to  QD  systems,  such  as  low  linewidth  enhancement  factors  and  reduced  values  of 
unsaturated  gain  and  absorption3’ 4,  if  understood  and  exploited  properly,  could  offer  considerable  performance  benefits 
over  bulk  and  quantum  well  (QW)-based  systems.  For  instance,  quantitative  knowledge  of  the  bias  conditions  under 
which  a  device  will  generate  stable  pulses  with  good  quality  is  highly  desirable  for  any  application  where  stable  pulses 
are  desired  under  variable  operating  conditions.  A  typical  example  of  such  an  application  is  the  possible  use  of  QDMLLs 
as  pulsed  sources  in  operating  environments  such  as  data  centers  (where  uncooled  QDMLLs  are  sought  as  sources  of 
stable  pulses)  that  require  components  to  be  able  to  tolerate  high  temperatures  accompanied  by  rapid  and  considerable 
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heat  generation  within  the  system.  Consequently,  devices  need  to  be  engineered  to  be  able  to  sustain  stable  output  pulses 
over  a  range  of  operating  conditions5, 6. 

In  order  to  design  and  engineer  devices  that  can  meet  the  stringent  restrictions  of  such  precision  applications,  robust 
modeling  techniques  are  highly  desirable.  And  while  a  considerable  amount  of  effort  has  gone  into  developing  detailed 
analytical  and  numerical  models  to  describe  mode  locking  in  semiconductor  lasers  (and  even  modified  specifically  for 
QDMLLs),  the  information  available  in  the  literature  about  device  parameter  values  unique  to  QD-based  structures  and 
how  they  can  be  accessed  and  exploited  so  as  to  optimize  device  perfonnance,  is  very  limited.  This  has  led  to  the 
associated  theoretical  and  experimental  work  following  two  distinct  paths,  with  the  former  focusing  on  developing 
analytical/numerical  models  of  varying  degrees  of  complexity  to  describe  specific  physical  processes/phenomena 
constituting  the  mode-locking  process,  while  the  latter  typically  relies  on  a  process  of  iterative  design  through  extensive 
testing  and  characterization.  On  the  theoretical  side,  as  more  physical  effects  and  processes  are  incorporated  into  the 
analysis,  the  model  grows  in  complexity  and  evolves  to  encompass  a  large  parameter  space.  Although  this  added 
complexity  should  enable  a  better  understanding  of  the  underlying  device  physics  and  operation,  the  lack  of  easily 
measurable  values,  coupled  with  the  number  of  required  parameters  often  limits  the  practical  utility  of  more  rigorous 
models. 

Previously,  we  used  a  Delay  Differential  Equation  (DDE)  model  for  passive  mode  locking  in  semiconductor  lasers7, 8' 9  to 
analyze  the  asymmetry  inherent  in  the  output  pulses  generated  by  QDMLLs10.  In  this  paper,  the  DDE  model  is  used  to 
study  the  impact  of  the  linewidth  enhancement  factor  in  the  gain  and  saturable  absorber  sections  on  the  steady  state 
device  output.  The  use  of  a  system  of  differential  equations  incorporating  a  delay  tern  as  opposed  to  a  system  of  partial 
differential  equations  enables  a  reduced  parameter  space  and  computational  simplicity11.  To  this  end,  physical  quantities 
measured  from  the  gain  and  loss  spectra  of  the  QD  substrate  material  are  transformed  into  model  parameters,  which  are 
then  used  as  initial  conditions  to  seed  the  model.  Relative  to  the  other  models  used  to  study  passive  mode  locking  in 
semiconductor  lasers,  the  DDE  model  has  a  parameter  space  of  just  9  parameters,  most  of  which  can  be  extracted  from 
measurable  physical  quantities,  leaving  only  a  small  number  of  free  parameters12, 13,  14.  This  enables  seeding  the  model 
with  realistic  input  parameters  for  each  operating  condition.  The  linewidth  enhancement  factors  in  the  gain  and  absorber 
sections,  which  are  not  available  from  the  gain  and  loss  data,  are  then  left  as  free  parameters  that  can  be  varied  to  study 
changes  in  device  output.  Getting  a  sense  of  the  linewidth  enhancement  factor  of  a  device  operating  under  a  certain  set  of 
conditions  is  important,  because  it  is  indicative  of  the  frequency  chirp  in  the  system,  and  QDMLLs  with  low  chirp  are 
highly  desirable  for  telecommunications  applications15.  Also,  numerical  analysis  has  shown  that  in  the  case  of  a  QDMLL 
operating  under  optical  feedback,  the  role  of  the  linewidth  enhancement  factor  in  the  gain  and  absorber  sections  is  a  key 
determining  factor  for  the  mode-locking  dynamics16. 

However,  experimental  measurement  techniques  are  known  to  yield  widely  varying  results,  depending  on  the  technique 
used  and  the  biasing  conditions17.  Consequently,  experiments  have  reported  values  ranging  from  very  small18  to  very 
large19.  For  this  reason,  it  is  highly  desirable  to  have  a  numerical  model  that  can  predict  an  estimated  range  of  values  that 
the  linewidth  enhancement  factor  can  take,  under  a  certain  set  of  operating  conditions. 

In  addition  to  giving  insight  into  the  range  of  permissible  values  that  the  linewidth  enhancement  factor  can  take,  the 
simulations  are  shown  to  predict  key  experimental  trends,  such  as  pulse  trimming  for  an  applied  bias  voltage.  The 
primary  objective  of  this  work  is  to  show  that  the  DDE  model,  used  in  its  original  form  without  any  modifications  and 
seeded  with  parameter  values  extracted  from  actual  gain  and  loss  measurements  on  the  constituent  QD  material,  provides 
a  highly  useful  guide  to  predict  and  analyze  device  perfonnance. 

2.  THEORETICAL  BACKGROUND 

The  Delay  Differential  Equation  (DDE)  model  due  to  Vladimirov  and  Turaev  is  derived  from  the  standard  coupled 
Partial  Differential  Equation  (PDE)  formalism,  which  describes  the  interaction  of  a  slowly  varying  complex  optical  field 
with  the  carrier  densities  in  the  gain  and  absorber  sections  of  a  generic  semiconductor  laser8, 9: 

dr 
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dr 

(3) 

In  equations  ( 1 ) — (3)  above,  ^incorporates  the  effect  of  spectral  filtering,  Q^and  cc  are  the  linewidth  enhancement 
factors  in  the  gain  and  absorber  media  respectively,  and  K accounts  for  linear  cavity  losses.  The  delay  parameter 
incorporates  history,  which  appears  in  the  first  equation  of  the  model  through  the  evaluation  of  optical  field  A ,  saturable 
gain  G ,  and  saturable  loss  Qat  the  times  X  —  T  ,  where  T  denotes  the  normalized  round  trip  time.  Parameter 

F  =  (ra6s  j Xgain  J  is  the  ratio  between  the  absorber  and  gain  relaxation  times,  and  S  is  the  saturation  parameter.  Finally, 
g0  and  q0 are  the  unsaturated  gain  and  absorption,  respectively. 

A  recent  experiment6  has  shown  that  the  segmented  contact  method20, 21  can  be  used  to  obtain  modal  gain  and  total  loss 
spectra  as  a  function  of  absorber  bias  voltage  and  temperature  for  the  QD  material  comprising  the  device.  These  curves 
were  used  to  extract  the  modal  gain  (gmod)  and  unsaturated  absorption  (ao)  as  a  function  of  current  density  at  each 
absorber  bias  condition.  In  the  following,  a  systematic  method  to  transform  physically  measured  parameters  from  the 
gain  and  loss  spectra  into  seed  conditions  for  the  DDE  model  is  discussed. 

Now,  from  equation  (2)  in  [6],  the  threshold  condition  for  lasing  can  be  expressed  as  follows, 


(gmod  (J)~ai  )Lg  -  ( a0  +  )La  =  «,  (4  +  4  ) 


(4) 


where,  am  and  OLi  represent,  respectively,  mirror  losses  and  internal  losses,  and  L,t  and  La  represent,  respectively,  the 


gain  and  absorber  section  lengths,  with: 


a~  =|  —  ]  In 


'  1  ^ 
v/44 


(5) 


where,  [L„  +  4  )  —  L  and  R{  ,  R 0  denote  mirror  reflectivities. 

In  the  DDE  model,  G  and  Q  describe  the  saturable  gain  and  loss  introduced  by  the  gain  and  absorber  sections, 


respectively,  and  K  <  1  describes  the  total  roundtrip  non-resonant  linear  intensity  losses.  Thus,  the  threshold  condition 
for  lasing  is  given  by9: 


Ke(G'Q)  =l=>G  =  g-ln(jc) 


Comparing  (3)  and  (4),  we  find: 


G 


threshold 


(gmod(*/)-«i)4 


(6) 


Qthreshold  ~  ( U0  +  01  i  )4 


(7) 


(8) 

k  =  ^R,R2 

(9) 

A  key  feature  of  the  DDE  model  is  its  dimensionless  formalism,  whereby  physically  measured  laser  parameters  must  be 
scaled  appropriately  and  made  dimensionless,  before  being  used  as  inputs  to  the  model.  Thus,  in  order  to  enable  a  direct 
conversion  between  measurable  static  laser  quantities  and  the  definitions  appearing  in  [9],  transformation  relations  to 
relate  the  saturable  gain  and  absorption  (G  and  Q,  respectively)  to  the  corresponding  unsaturated  parameters  (  goand  9o  > 
respectively)  have  been  derived  to  yield: 


g0  =  rG 


(10) 
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%  =  (2A) 


(ii) 


where, 


„  =  gf*  _  [^mod^Va/LuH. 


ggrg  [d8™AJ)l8A 


(12) 


Further,  using  the  dynamical  stability  analysis  presented  in  Fig.  5  of  [9],  it  is  deduced  that  a  suitable  constraint  to 
achieve  stable,  fundamental  mode-locking  in  terms  of  g0  and  q0,  yields  an  expression  for  the  carrier  relaxation  ratio, 
r,  which  depends  solely  on  measurable  static  device  parameters  defined  above  in  equations  (7),  (8)  and  (12): 


(13) 


Equations  (7)  -  (13)  provide  a  set  of  initial  conditions  that  can  be  used  to  seed  the  simulation. 


3.  SIMULATION  RESULTS 


In  this  section,  the  simulation  results  for  an  8-stack  InAs  dots-in-a-well  (DWELL)  passively  mode  locked  laser,  with 
a  1  mm  absorber  section  and  7  mm  gain  section  are  presented.  The  details  of  device  structure,  experimental  data  and 
results  can  be  found  in  [6],  The  approach  involves  extracting  parameters  appearing  on  the  right  sides  of  equations 
(7),  (8),  (12)  and  (13)  from  quantities  measured  as  a  function  of  gain-section  current  density,  wavelength,  absorber 
reverse  bias  and  temperature  using  the  segmented  contact  method20' 21 . 

This  method  of  parameter  extraction  provides  an  extremely  convenient  technique  to  determine  model  parameters  for 
input  into  the  DDE  model.  Such  a  situation  may  be  envisioned  for  a  case  where  device  performance  needs  to  be 
simulated  over  a  range  of  temperatures  and  operating  conditions.  In  such  cases,  the  benefits  of  being  able  to 
accurately  simulate  the  output  characteristics  of  a  proposed  device,  with  no  more  than  a  relatively  simple  set  of  gain 
and  loss  measurements  on  a  test  structure  (and  parameters  extracted  from  these),  and  using  the  results  to  guide 
device  engineering  can  be  clearly  seen  to  be  invaluable.  Conversely,  the  simulation  results  can  offer  invaluable 
insight  into  the  range  of  values  a  given  parameter  can  take,  in  order  for  device  operation  to  be  stable.  For  instance,  if 
the  values  extracted  for  the  unsaturated  gain  and  absorption  are  raised  by  a  factor  of  ten  each  (while  all  other 
parameter  values  are  retained),  the  simulation  yields  unstable  pulses.  In  this  section,  simulation  results  obtained  by 
seeding  the  model  with  parameters  extracted  from  physically  measured  quantities  are  presented  and  discussed. 

A  close  examination  of  the  system  of  DDEs  given  in  equations  (l)-(3)  reveals  the  parameter  space  involved.  Thus, 
in  order  to  model  a  QDMLL  realistically,  it  is  imperative  to  constrain  as  many  of  the  parameters  as  possible  to 
values  typical  of  a  QD  device.  Since  the  onset  of  lasing  occurs  at  threshold,  the  first  step  is  to  use  values  measured  at 
threshold  in  equations  (7)  and  (8)  to  calculate  G(0)  and  Q(0),  and  then  use  the  values  of  G(0)  and  Q(0)  in  equations 
(10)  and  (11)  to  calculate  the  corresponding  unsaturated  parameters,  g0  and  q0  at  threshold.  Further,  parameters  k 
and  s  can  be  directly  found  from  equations  (9)  and  (12)  respectively,  while  the  carrier  relaxation  ratio  (T)  is 
extracted  from  equation  (13).  The  delay  parameter  T  is  determined  for  each  set  of  operating  conditions  by  scaling 
the  cavity  round  trip  time  to  the  corresponding  absorber  recovery  time,  which  is  chosen  as  the  characteristic 
timescale.  Of  the  remaining  three  model  parameters,  the  spectral  filtering  coefficient  (y)  was  judiciously  chosen  and 
held  constant,  while  the  values  of  linewidth  enhancement  factor  in  the  gain  and  absorber  sections  were  varied  in 
tandem.  Thus,  the  main  assumptions  made  to  simplify  the  analysis  were:  1 .)  the  spectral  filtering  was  assumed  to  be 
invariant  with  respect  to  the  gain  and  linewidth  enhancement  factors,  and  2.)  the  linewidth  enhancement  factors  in 
the  gain  and  absorber  sections  were  kept  equal. 

The  results  shown  in  the  following  were  obtained  as  follows.  The  parameter  values  measured/extracted  as  discussed 
above  were  used  as  initial  conditions  for  the  simulations.  The  system  of  equations  defining  the  model  (equations  (1), 
(2)  and  (3))  was  then  integrated  over  several  round  trip  times  to  allow  transients  to  settle.  The  primary  objective  was 
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to  constrain  as  many  of  the  model  parameters  as  possible  to  measured  values,  so  as  to  get  a  sense  of  the  range  of 
values  that  the  non-measured  parameters  (namely,  the  linewidth  enhancement  factors  in  the  gain  and  absorber 
sections)  can  take  under  the  operating  conditions  under  study.  Stable  device  operation  is  simulated  as  stable, 
fundamental  mode  locking  (identical  pulses,  once  every  round  trip),  so  that  any  variation  in  pulse  amplitude, 
duration  or  structure,  or  multiple  pulses  in  a  round  trip  are  representative  of  a  breakdown  of  the  device  output. 

Table  I  lists  a  sample  set  of  parameter  values,  extracted  from  gain  and  loss  spectra  at  20  C  and  no  applied  absorber 
bias  at  threshold  (current  density  =  357  A/cm2).  The  delay  parameter  T  has  been  scaled  using  a  characteristic 
timescale  of  62  ps,  which  was  inferred  from  measurements  performed  on  a  QD  saturable  absorber  with  a  similar 
epitaxial  structure22,  using  an  equivalent  value  of  active  region  field  strength. 


Table  1:  Simulation  Parameters  at  20  C,  0  V  Absorber  Bias  (qg  =  a0  =  2.5) 


T=  3.226 

G(0)  =3.33 

g(0)=2.33 

«g=2.5 

a,  =  2.5 

S  =  2.676 

r=o.i3 

y=  13.3 

xr  =  0.55 

The  values  of  OC„  and  a  shown  are  for  one  particular  case.  They  have  been  chosen  equal  based  on  the  analysis 

presented  in  [9],  where  it  was  numerically  shown  that  the  best  pulse  quality  is  achieved  for  the  case  of  equal  values 
of  linewidth  enhancement  factor  in  the  gain  and  absorber  sections. 


Fig.  lashowsthe  steady  state  output  of  the  system  of  equations  in  (l)-(3),  for  the  case  OC  g=  aq  =  2.5. 

In  the  absence  of  an  applied  bias  on  the  absorber  section,  there  is  no  sweep-out  mechanism  as  a  result  of  the 
associated  electric  field  to  expedite  the  carrier  relaxation  rate  in  the  absorber  section.  Thus,  the  absorber  takes  a 
longer  time  to  recover,  resulting  in  sub-optimal  pulse  trimming.  Thus,  the  simulation  results  predict  broad  pulses, 
with  a  width  ~  23  ps  (Fig.  lb). 
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Fig.  la:  Steady  state  output  for  0  V  absorber  bias,  ag  =  aq  =  2.5:  G(x),  Q(x),  (G(x)-Q(x))  and  |A(x)|2. 

STEADY  STATE  FIELD  INTENSITY 
(T  =  20  C,  J  =  Jth  =  357  A/cm2,  Vabs  =  0  V,  alpha  =  alpha  =  2.5) 
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Further,  the  mode-locked  pulses  are  predicted  to  become  unstable  for  ag  =  a  =  2.9  (Fig.  lc).  This  is  evidenced  by 


the  unequal  pulse  intensities,  along  with  the  appearance  of  a  pronounced  side  structure  on  the  leading  edge  of  each 
pulse. 


STEADY  STATE  FIELD  INTENSITY 
(T  =  20  C,  J  =  Jth  =  357  A/cm2,  Vabs  =  0  V,  alphag  =  alphag  =  2.9) 


It  may  be  noted  that  mild  instabilities  were  obtained  for  the  cases  where  the  a-parameters  were  chosen  to  be  greater 
than  2.5.  The  case  CCg  =  CCq  =  2.9  is  shown  here  to  elucidate  the  nature  of  typical  instabilities. 

Thus,  the  simulation  results  show  that  for  the  combination  of  parameters  extracted  at  20  C  with  no  applied  reverse 
bias,  the  values  of  linewidth  enhancement  factor  in  the  gain  and  absorber  sections  must  be  less  than  or  equal  to  2.5, 
in  order  to  ensure  stable,  fundamental  mode-locking. 

Next,  the  effect  of  an  applied  (moderate)  bias  voltage  is  considered,  with  Vabs  =  -3V.  Again,  in  this  case,  all  other 
parameters  are  constrained  to  extracted  values,  and  the  effect  of  the  linewidth  enhancement  factor  is  analyzed  by 
simulation.  Table  II  lists  a  sample  set  of  parameter  values,  extracted  from  gain  and  loss  spectra  at  20  C  and  an 
applied  absorber  bias  of  -3V  at  threshold  (current  density  =  462  A/crn2).  Again,  using  the  measurements  reported  in 
[22],  the  delay  parameter  T  was  scaled  to  a  characteristic  timescale  of  40  ps. 


Table  II:  Simulation  Parameters  at  20  C,  -  3  V  Absorber  Bias  (qe  —  a0  =  1.8) 


r=5.o 

G(o)  =4.18 

g(0)  =  3.20 

«g=L8 

aq  =  1  -8 

S  =  4.65 

r=o.o8 

y=  13.3 

K=  0.55 

Fig.  2a  shows  the  simulation  output  for  the  case  OL  =  (X  -  1.8.  Compared  to  the  previous  case  with  no  applied 

bias,  the  net  gain  profile  in  the  third  panel  shows  a  sharper  profile,  suggesting  a  shorter  interval  of  high  net  gain  - 
this  is  seen  as  considerably  better  trimmed  pulses  in  the  fourth  panel,  with  a  pulse  width  ~  10  ps  (Fig.  2b). 


STEADY  STATE  OUTPUT 

(T  =  20  C.  J  =  Jlh  =  462  A/cm7,  Vabs  =  -3  V,  alpba^  =  alplw,  =  1.8) 
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Fig.  2a:  Steady  state  output  for  -3  V  absorber  bias,  ag  =  aq  =  1.8:  G(x),  Q(x),  (G(x)-Q(x))  and  |A(x)|2. 
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Fig.  2b:  Steady  state  pulse  intensity  profile  for  -3  V  absorber  bias,  ag  =  aq=  1.8. 

In  this  case,  the  simulation  results  predict  that  instabilities  begin  to  appear  in  the  mode- locked  pulses  for  CC g  =  CCq  - 
1.9  (Fig.  2c),  as  seen  from  the  variable  pulse  intensities. 

STEADY  STATE  FIELD  INTENSITY 
(T  =  20  C,  J  =  Jth  =  462  A/cm2.  Vabs  =  -3V,  alpha  =  alpha  =  1 .9) 


Finally,  we  consider  the  effect  of  a  strong  applied  bias  of  -5  V  on  the  absorber  section.  Table  III  lists  a  sample  set  of 
parameter  values,  extracted  from  gain  and  loss  spectra  at  20  C  and  an  applied  absorber  bias  of  -5V  at  threshold 
(current  density  =  537  A/cm2).  Again,  using  the  measurements  reported  in  [22],  the  delay  parameter  T  was  scaled  to 
the  characteristic  timescale,  chosen  to  correspond  to  an  absorber  recovery  time  of  30  ps. 


T=  6.667 

G(0)  =4.55 

Q{  0)=3.62 
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T  =  0.0576 

y=  13.3 

*-  =  0.55 

Fig.  3a  shows  the  simulation  output  for  the  case  Clu  =  OCq  =  0.8. 

In  this  case,  we  find,  comparing  the  first  and  second  panels  in  Fig.  3a  to  the  corresponding  panels  in  Fig.  2a,  that 
while  the  saturable  gain  profile  in  both  cases  follows  a  similar  trend,  the  saturable  loss  saturates  at  a  lower  value  for 
a  higher  bias.  The  lower  losses  lead  to  a  slightly  higher  net  gain  window  than  the  -3  V  case,  so  that  pulses  are 
trimmed  to  ~  8.75  ps  (Fig.  3b). 
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STEADY  STATE  OUTPUT 

(T  =  20  C.  J  =  Jth  =  537  Alcm7,  Vabs  =  -5  V.  alphag  =  alphaq  -  0.8) 
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Fig.  3a:  Steady  state  output  for  -5  V  absorber  bias,  ag  =  aq  =  0.8:  G(x),  Q(x),  (G(x)-Q(x))  and  |A(x)|2. 

STEADY  STATE  FIELD  INTENSITY 
(T  »  20  C,  J  =  Jth  *  537  AJcvn2,  Vabs  =  -5  V,  alpha  =  alpha  =  0.8) 


Fig.  3b:  Steady  state  pulse  intensity  profile  for  -5  V  absorber  bias,  tig  =  aq  =  0.8. 


Interestingly,  the  onset  of  instabilities  in  this  case  appears  at  a  =  CLq  =  0.9  as  an  abrupt  switch  from  fundamental 

to  second-order  mode-locking  (Fig.  3c).  This  change  of  operating  regime  is  characterized  by  the  sudden  appearance 
of  two  pulses  every  round  trip.  Physically,  this  represents  a  breakdown  of  stable,  fundamental  mode-locking. 


STEADY  STATE  FIELD  INTENSITY 
(T  *  20  C,  J  a  Jth  «  537  A/cm3,  Vabs  *  -5V,  alpha  =  alpha  =  0.9) 


Fig.  3c:  Steady  state  field  intensity  |A(r)|2  for  -5  V  absorber  bias,  ag  =  aq  =  0.9. 


A  comparison  between  Fig.  la,  2a  and  3a  shows  that  the  impact  of  increasing  absorber  bias  voltage  is  a  slight 
sharpening  of  the  saturable  gain  and  loss  profiles,  so  that  a  slightly  higher  peak  gain  coincides  with  a  slightly  lower 
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loss  in  the  latter  case,  resulting  in  larger  net  gain  in  a  shorter  time  window.  This  is  predicted  to  lead  to  narrower 
pulses. 

Further,  the  application  of  a  bias  voltage  not  only  trims  the  output  pulses  considerably,  but  also  results  in  a 
decreased  range  of  possible  values  of  the  linewidth  enhancement  factor  in  order  to  maintain  stable  mode-locking. 
Thus,  the  range  of  possible  values  of  GQand  a  extends  to  2.5  when  no  reverse  bias  is  applied,  but  is  restricted  to 
1 .8  or  less  when  a  moderate  bias  of -3  V  is  applied,  and  0.8  or  less  when  a  strong  bias  of  -5  V  is  applied. 

4.  CONCLUSIONS 

A  Delay  Differential  Equation-based  approach  has  been  used  to  model  a  two-section,  passively  mode-locked 
quantum  dot  laser,  where  most  of  the  model  parameters  have  been  extracted  from  gain  and  loss  measurements 
performed  on  the  quantum  dot  substrate.  The  objective  was  to  get  a  sense  of  the  linewidth  enhancement  factor  in  the 
gain  and  absorber  sections  distinctive  to  the  particular  QD  device  operating  under  a  certain  set  of  conditions,  in 
order  to  ensure  stable,  fundamental  mode-locking.  To  this  end,  the  parameters  that  were  not  obtainable  from  the 
gain  and  loss  measurements,  namely,  the  linewidth  enhancement  factors  in  the  gain  and  absorber  sections,  were 
varied  and  their  impact  on  the  device  output  was  studied  at  a  constant  temperature  of  20  C,  for  the  cases  of  no 
applied  bias,  moderate  bias  (-  3V)  and  strong  bias  (-  5V)  on  the  absorber  section.  The  key  simplifying  assumptions 
were  to  keep  the  spectral  filtering  coefficient  invariant  with  respect  to  the  gain  and  linewidth  enhancement  factors, 
and  to  consider  equal  values  of  linewidth  enhancement  factors  in  the  gain  and  absorber  sections.  The  simulation 
results  not  only  predicted  the  experimentally-observed  trend  of  pulse  narrowing  with  applied  bias,  but  also  that 
stable  mode  locking  is  restricted  to  a  narrow  range  of  possible  values  of  the  linewidth  enhancement  factors,  before 
either  the  pulses  become  unstable,  or  the  device  switches  to  higher-order  mode  locking.  Further,  the  results  showed 
that  this  range  is  significantly  reduced  with  the  application  of  a  moderate  bias  (and  reduced  even  further  for  a  strong 
bias),  which  is  important  for  telecom  and  datacom  applications,  as  well  as  for  military  RF  photonic  applications 
such  as  low  noise  clocks. 
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